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A  turbofan  engine  propulsion  system  in  which  large  amounts  of  water 
are  injected  into  the  fan  discharge  duct  is  investigated  with  the  goal 
of  increasing  both  the  thrust  and  propulsive  efficiency  while  retaining 
the  light-weight  qualities  of  a  standard  turbofan  engine.   A  parametric 
computer  analysis  is  used  to  examine  the  effect  of  several  variables, 
including  water-to-gas  generator  air  ratio.,  water  injection  velocity, 
fan  duct  pressure  loss,  and  fan  duct  thermal  and  dynamic  nonequilibrium, 
upon  thrust  and  propulsive  efficiency.   In  addition,  the  design  param- 
eters of  fan  pressure  ratio  and  fan  bypass  ratio  are  examined  for  their 
optimum  values,  and  optimum  operating  combinations  of  water-to-gas  ratio 
and  water  injection  velocity  are  determined. 

A  test  apparatus  is  developed  for  the  direct  measurement  of  wall 
friction  force  in  two-phase  flows.   A  computer  program  is  presented  to 
reduce  experimental  data  and  compare  with  pressure  drop  predicted  by 
two  empirical  correlations. 
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I.   COMPUTER  OPTIMIZATION  OF  WATER-AUGMENTED  TURBOFAN  CONCEPT 

A.   SOURCE  AND  BACKGROUND 

Transportation  by  water-borne  craft  is  one  of  the  oldest  and  most 
economical  means  known  to  man.   Although  many  improvements  have  been 
made  in  payload  carrying  ability,  navigational  ability  and  propulsion 
systems  over  the  centuries,  the  degree  of  improvement  in  the  speed  of 
water  craft  has  been  quite  small.   It  has  become  desirable,  from  both 
a  military  and  civilian  point  of  view,  to  find*  new  light-weight  pro- 
pulsion systems  for  water  craft  which  provide  significant  increases  in 
speed  capability  while  retaining  the  load  carrying  ability  and  effi- 
ciency of  previous  sea  transportation  systems. 

The  feasibility  of  producing  a  ship  or  boat  capable  of  sustained 
high  speeds  with  relatively  large  payloads  and  reasonable  efficiency 
depends   largely  upon  two  factors:   drag  reduction  and  the  availabil- 
ity of  a  light-weight  efficient  propulsion  system.   Drag  reduction 
can  be  achieved  by  use  of  hydrofoils,  the  captured  air  bubble  concept 
or  recent  boundary  layer  control  schemes  including  blowing,  suction 
and  polymer  injection.   The  development  of  a  suitable  propulsion  sys- 
tem, however,  has  not  be  quite  so  active. 

The  turbofan  provides  relatively  high  thrust  and  very  low  system 
weight  compared  to  more  conventional  marine  propulsion  systems.   How- 
ever, in  the  speed  range  contemplated  for  marine  craft,  100  knots  and 
below,  the  propulsive  efficiency  of  the  turbofan  is  too  low  to  allow 
serious  consideration. 

Early  methods  of  increasing  the  thrust  and  propulsive  efficiency 
in  turbojet  and  turbofan  engines  centered  on  injection  of  water  into 
the  compressor  or  exhaust  section.   This  method  of  water  injection 
has  been  used  primarily  in  aircraft  applications  for  short  periods 
only  and  has  several  drawbacks  for  marine  applications,  among  which 
are  stiff  requirements  on  water  purity  due  to  machinery  corrosiveness 
and  energy  transfer  losses  caused  by  water  evaporation  due  to  mixing 
with  hot  gases.   A  more  recent  proposal  has  been  to  inject  large  quan- 
tities of  water  into  the  fan  discharge  duct  of  a  turbofan  engine,  thus 
overcoming  the  two  primary  drawbacks  mentioned  above  because  the  gases 
are  relatively  cool  and  the  water  contacts  no  moving  machinery.   Con- 
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sequently,  sea  water  can  be  used  in  large  quantities. 

Several  investigators  have  made  preliminary  studies  concerning  the 
feasibility  of  using  water-augmented  air-jets  for  water  craft  propul- 
sion.  Muench  and  Keith   IRef.  lj   conducted  a  preliminary  parametric 
study  of  possible  propulsive  efficiencies  and  concluded  that  reason- 
able efficiencies  can  be  realized.   Davison  and  Sadowski   Ref .  2J 
applied  their  analysis  to  a  particular  existing  turbofan  engine  design 
with  water  injected  into  the  fan  discharge  duct.   Their  analysis  showed 
that  although  efficiencies  were  slightly  lower  than  for  other  more  con- 
ventional marine  propulsion  systems,  the  extremely  low  system  weight 
more  than  compensated  for  the  efficiency  deficit.   Quandt   Ref.  3| 
has  assembled  some  recent  theoretical  and  experimental  techniques  con- 
cerning "the  gas-phase-continuous  two-phase  jet  system." 

Knudson   IRef.  4j   followed  the  work  of  Davison  and  Sadowski  in 
proposing  the  injection  of  large  amounts  of  sea  water  into  the  fan  dis- 
charge duct  of  a  turbofan  engine  to  achieve  increased  thrust  and  pro- 
pulsive efficiency.   The  general  arrangement  of  the  water-augmented 
turbofan  engine  is  presented  in  Fig,  1.   Knudson  was  the  first  to  con- 
sider the  effects  of  velocity  slip  ratios  and  temperature  ratios  be- 
tween water  droplets  and  air  in  the  fan  mixing  duct  and  nozzle.   The 
effects  of  variation  of  water-  to-gas  mass  flow  ratio  and  water  injection 
velocity  were  also  examined. 

The  present  analysis  undertakes  a  more  detailed  computer  study 
which  includes  the  pump  work  required  to  inject  the  water  at  velocities 
greater  than  craft,  velocity.   A  wide  range  of  design  and  operating  param- 
eters is  examined  for  the  effect  on  overall  thrust  and  propulsive  effi- 
ciency.  Some  minor  programing  errors  were  discovered  in  Ref.  4  and 
these  are  corrected  in  the  present  analysis. 

B.   WATER- AUGMENTED  TURBOFAN  ANALYSIS 

The  analysis  which  follows  is  based  on  that  given  by  Knudson  with 
appropriate  changes  and  corrections  for  the  present  use.   The  analysis 
differs  from  that  for  a  normal  "dry"  t.urbofan  in  the  fan  mixing  duct 
and  fan  nozzle  sections  due  to  the  injection  of  water.   The  equations 
are  developed  in  a  manner  similar  to  their  use  in  the  computer  program. 


12 


1 .  Assumptions 

Following  general  practice  in  the  analysis  of  turbofan  engines 
the  air  was  assumed  to  behave  as  a  perfect  gas  with  two  constant  values 
of  specific  heat  ratio,  d  ,  1.4  for  the  relatively  cool  sections  of  the 
system  and  1.33  for  the  high  temperature  sections  of  the  engine.   In  the 
burner  section  of  the  engine  the  arithmetical  mean  of  the  two  values  of 
specific  heat  ratio  was  used  since  the  air  transitioned  from  the  cool  to 
the  hot  state.   It  was  assumed  that  no  heat  was  transferred  except  in 
the  combustion  chamber  and  that  no  work  was  transferred  other  than  in 
the  compressor,  fan  and  turbine.   The  mixing  duct  area  was  assumed  con- 
stant and  it  was  also  assumed  that  only  converging  nozzles  were  used  in 
the  engine. 

Finally,  it  was  assumed  that  completely  dry  air  entered  the 
system.   The  external  drag  of  the  water  inlet  scoop  was  not  included 
in  the  analysis  since  this  factor  is  dependent  on  the  actual  scoop 
design  employed. 

2 .  Gas  Generator  Diffuser 

The  ambient  total  pressure  is  determined  from 

PTO  ■  K.U  *  ^Ml)1^  (i) 

and  the  diffuser  efficiency  is  defined  by 


'     Pro   ~  Pao 

Thus 


(2) 


Pr±    =  Pao    *   $,  (.Pro  '  Pao)  (3) 

The  total  temperature  remains  constant  through  the  diffuser  and  is 

The  total  enthalpy  is 

Hri  =     Cpc(Tti  '  Tao)   +    Hfi£pA  (5) 
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where  C   is  the  cool  temperature  specific  heat  and  H     is  the 

enthalpy  of  the  air  at  the  ambient  static  temperature,  T.^. 

J      AO 

3.   Compressor 

The  ideal  (isentropic)  temperature  at  the  compressor  exit  for 

a  specified  total  pressure  ratio,  P  „/P    is 

iTc-1 

'c  (6) 


Trzi   -  TT1(^/Prt)^ 


If  the  compressor  efficiency,  J7q     ,  is  defined  as  the  ratio  of  the 
isentropic  compressor  work  to  the  actual  compressor  work,  the  total 
enthalpy  at  the  compressor  exit  is 

u  CpcCTrzz- TTi)  u 

rtrz  -  ^  ^  riTi  (7) 

and  the  corresponding  total  temperature  is 

4 .   Burner 

The  specific  heat  in  the  burner  section  is  taken  as 

/"»         Cpc  +    Cph 

^PB    =       ^ (9) 

The  burner  exit  total  enthalpy  is 

Hr3   =  CpB(TT3-TTZ)  t  Htz  <io) 

where  the  total  temperature  at  the  burner  outlet,  T     is  fixed  by  the 
maximum  allowable  turbine  inlet  total  temperature.   The  energy  equation 
through  the  burner  can  be  written 


rhc  Htz  *  »%  (HfnZ  ^bHVf)    =   (^c^^fj  ^3 


(ii) 


where  iiu  is  the  compressor  mass  flow  rate,  mr      is  the  fuel  mass  flow 
C  fu 

rate,  H    is  the  fuel  enthalpy  prior  to  entering  the  burner,  r) g    is 
the  burner  efficiency,  HVF  is  the  heating  value  of  the  fuel  and  H   is 
the  total  enthalpy  at  the  burner  outlet.   If  it  is  assumed  that  the 
fuel  injection  temperature  is  equal  to  the  air  temperature  at  the 
injection  point  and  that  the  energy  needed  to  vaporize  the  fuel  is 
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included  by  choosing  the  lower  heating  value  of  the  fuel,  then  the 
fuel-air  ratio  can  be  represented  by 

FAR  -     *S     =    Crh(rr}-Tr2) 

5    Turbine 

The  total  work  output  of  the  turbine  must  equal  the  total  work 
input  to  the  other  components  of  the  engine . 

where  m  is  the  turbine  mass  flow  rate,  itl,  is  the  fan  mass  flow  rate, 

m„  is  the  water  mass  flow  r; 
w 

work  per  unit  mass  given  by 


mr,  is  the  water  mass  flow  rate  and  W„  is  the  water  infection  pump 
W  P 


2  ,„-     ..  ,Z 


Vjr    -C*«Vo) 


Wp  =    -  2^;—  en) 


where  VrT_,  and  V~  are  the  water  iniection  velocity  and  the  craft 
W7      0 

velocity,  respectively,  7/^   is  the  water  intake  system  velocity  recov- 
ery factor  and  *T)p     is  the  pump  efficiency.   The  total  temperature  at 
the  turbine  outlet  is 

<T4      =  lT3   —   ^ (14) 

If  the  turbine  efficiency  is  defined  as  the  ratio  of  actual  work  output 
to  isentropic  work  output,  then 

7~  _     J"  Ht3  -    Mt4 

IT4T    -     ir*  —     ^ 


and 


T4t  =  'T3 ?,<:,<■  (15) 


Pt*     =    Pr.r    -    Pr'(r"Z/rTt)^ri 


T.)    ■    n-n    -    tr3  /T  (16) 


6 .   Gas  Generator  Nozzle 

The  nozzle  must  discharge  to  atmospheric  pressure  (unless 
choked)  and  thus  the  ideal  static  temperature  at  the  nozzle  exit  is 


75«   •  r„  (%•«)    *h  <17> 
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Under  assumption  of  constant,  specific  heats,  the  nozzle  efficiency  can 
be  represented  by 

Tr"  "  T"S  (18) 


7,  =    »    p 


Thus,  for  a  given  nozzle  efficiency,  the  static  temperature  at  the 

nozzle  exit,  T.,_,  can  be  found  and  the  nozzle  exit  velocity  can  be 

'   A5  J 

determined  from 


to-  =   i2Cph(rrs-7js)' 


(19) 


The  exit  mach  number  is  given  by 


**  Vas        _  Mas 

Ws  =   — =       ,  ==r  (20) 

If  the  exit  mach  number  exceeds  unity,  the  computations  are  adjusted 
to  produce  sonic  velocity  at  the  exit  with  the  appropriate  exit 
pressure . 

7 .   Fan  Diffuser  and  Fan 

The  fan  diffuser  analysis  is  identical  to  that  of  the  gas  gen- 
erator diffuser.   Similarly,  the  fan  analysis  corresponds  to  the  com- 
pressor analysis. 

8 .   Fan  Mixing  Duct  Water  Injection  Analysis 

The  analysis  of  the  water  injection  process  in  the  fan  mixing 
duct  considers  the  possibility  of  variation  in  the  thermal  and  dynamic 
equilibrium  in  the  fan  mixing  duct  and  the  fan  nozzle,  respectively, 
and  static  pressure  changes  between  the  water  injection  plane  (fan 
mixing  duct  entrance)  and  the  fan  mixing  duct  exit  (stations  7  and  8, 
respectively,  in  Fig.  1). 

The  following  terms  are  used  throughout  the  analysis  and  are 
defined  here  for  convenience: 

Bypass  Ratio  (BR)  =  ""tf  An  £ 

Water  to  Gas  Generator  Air  Ratio  (WGR)   =  7nv/ffi^> 
Mixture  Ratio  (MR)  =  Wv /Trip 
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The   mean   velocity   of    the    two-phase    flow   is    defined   by 

and    the   mean  velocity    is    thus 

V     ~  liMR 

The  momentum  equation  in  the  mixing  duct  is 


(~mF  i  >h„)dv  =  -  AJP 


(21) 

If  the  area  occupied  by  the  air  is  m   / p  .V .  and  the  area  occupied  by 
the  water  is  A_-/a,V_,,  the  total  flow  area  can  be  written 

**^%  +  -gk)  <22> 

Combining  equations  (21)  and  (22)  and  assuming  a  constant  area  mixing 
duct  the  resulting  equation  can  be  integrated  yielding 


vt  ■  ^(^-^(t^JM 


MR 


flnV* 


/»T 


(°«-iV* 


wt 


(23) 


In  order  to  conduct  an  energy  analysis  through  the  mixing  duct 

it  is  necessary  to  consider  the  effects  of  vaporization  of  water. 

The  values  Hr,  H,_  ,  Sr  and  Sr      were  determined  for  the  computer  using 
f    fg    f       fg 

least  squares  cubic  approximations  of  the  data  in  the  steam  tables 

[Ref.  5j   for  the  temperature  range  510-660  degrees  Rankine .   The 

equations  used  are  shown  at  the  beginning  of  the  computer  program  in 

Appendix  II.   The  values  of  H  and  S   are  found  from  the  equations 

g      g 

H  =  H^  +  H 
g    f    fg 

s  =  sr  +  S^ 
g       f       fg 

The  partial  pressure  of  water  vapor,  Ppv,  in  atmospheres  is 
determined  from  equation   (12)  in  Ref.  5. 

A  ■ 


,09*>(ppv)s  l°3«>(pc)-  T'    A  ti  +  Dx "J 


(24) 


The  specific  humidity,  X,  of  the  mixture  is  defined  as  the  ratio  of  the 
mass  of  the  water  vapor  to  the  mass  of  the  dry  air  in  the  duct.   Using 
the  perfect  gas  relations  for  air  and  water  vapor  it  can  be  shown  that 


X8  -  (- 


PV8 


fAWy 


Taa 


Pa8-Ppv8   I  KM*       \  T*B 


(25) 
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where  MW  is  the  molecular  weight  and  subscripts  V  and  A  signify  water 
vapor  and  air,  respectively.   The  vapor  temperature  is  assumed  equal 
to  the  air  temperature. 

The  energy  equation  through  the  mixing  duct  is 

<26) 

The  analysis  thus  consisted  of  the  equations  (23),  (25)  and  (26)  with 

six  unknowns:   VrTO,  VAO,  Trio,  TAO,  P  _  and  X0  .   A  solution  is  found 
Wo    Ao    Wo    Ao    Ao       o 

by  specifying  values  for  three  ratios  of  unknowns: 

TR8WA   =   TW8/TA8  (27) 

VR8WA  =  Vw8/VAg  (28) 

PR87   =   PA8/PA7  (29) 

By  varying  these  ratios,  the  effects  of  departure  from  thermal  or 
dynamic  equilibrium  at  the  mixing  duct  exit  and  static  pressure  change 
through  the  duct  can  readily  be  examined.   Although  six  equations  with 
six  unknowns  exist,  an  explicit  solution  is  not  possible  because  equa- 
tions (25)  and  (26)  contain  involved  functions  of  Tr70  through  the  en- 

Wo 

thalpy  and  partial  pressure  terms.   However,  it  is  possible  to  deter- 
mine VAO  and  VrTO  using  equations  (23),  (28)  and  (29).   The  remaining 

Ao       Wo 

three  equations,  (25),  (26)  and  (27),  are  then  manipulated  to  obtain 
a  function  of  T.„  only  (called  FUN).   A  Newton-Raphson  iterative  solu- 
tion is  then  used  in  the  form 

where  subscript   j  indicates  the  j    approximation,  j-1  indicates  the 
(j-1)    approximation,  etc.,  and  DFUN  is  the  derivative  of  FUN  with 


res 


pect  to  TAO,  both  functions  being  evaluated  at  TAO .  . .   The  correct 
K        A8'  A8j- 1 


value  of  T.n    is  obtained  when  the  absolute  value  of  FUN /DFUN  becomes 
A8 

zero.   It  is  then  a  simple  matter,  knowing  T"   ,  to  determine  X  and 

T   from  equations  (25)  and  (27). 
Wo 


9 .   Fan  Nozzle  Analysis 

Following  reasoning  similar  to  that  in  the  preceding  section 
the  specific  humidity  at  the  fan  nozzle  exit  is 

The  constant  area  assumption  made  to  allow  solution  of  the 
momentum  equation  for  the  fan  mixing  duct  analysis  obviously  does  not 
apply  through  the  nozzle.   The  approach  taken  (as  in  Ref.  2)  is  that 
the  net  change  in  entropy  through  the  nozzle  is  zero. 

(S9  "  VaIR  +  <S9  "  VWATER   "   °  (32) 

Equation  (32)  when  expanded  gives 

t'[(*w-  Ws«  * »,  y,  fy]  -  o        <33) 

The  exit  pressure,  PAQ,  must  equal  atmospheric  pressure,  leaving  two 
equations,  (31)  and  (33),  with  unknowns  T    T   and  X    By  specifying 

TR9WA   =  t"Vta9  (34) 

and  solving  equations  (31)  and  (33)  for  a  function  of  T  _,  a  Newton- 
Ay 

Raphson  iteration  again  provides  the  solution  for  T  _,  T   and  X_ 
The  energy  equation  is 

=  M***^*(^-X*$fa*!()*»M6k*l£)  (35) 


The  fan  nozzle  exit  velocity,  V..,  is  determined  by  specifying 

VR9WA  =  VW9/VA9  (36) 

and  solving  equation  (35)  for  V._  as  shown  on  the  following  page. 

Ay 
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V 


CPc(Tflg-TAg)  4  XaHfjs-  X9Mf39    f  MR(^6-Hfq) 


fc[l4  X9*CMR-X9XVR9^)^ 
+      frill*  *»»  IflR-  X»)(  V/?5  W/Q*j 

[i  +x9+Uf?-x9XVR.9^)ij 


(37) 


10.      Specific   Thrust   and   Propulsive    Efficiency 

The    thrust   of    the   water-augmented    turbofan   is   determined    from 
the  momentum  flux  equation 

TH    =    -mTVfl5  -  Wc\/<,   +  (inf  +  X9TnF)VA^   -   yr\f  ]/0 


+■  (-mw7-X9^F)VW9-  ^W7V„,0  +  A5(pAS-PAo) 


38) 


where  the  water  velocity  entering  the  inlet  scoop,  V  _,  is  assumed 
equal  to  the  craft  velocity,  V_ .   The  specific  thrust  is  the  thrust 
divided  by  total  air  mass  flow  through  the  engine. 

-B1?lV.)  +  tW6R-X9BR)Vw9-WSR(»/w.>-^^-(V?|  39> 

The  propulsive  efficiency  is  a  measure  of  how  well  the  system  converts 
the  kinetic  energy  change  of  the  working  medium  into  thrust  power  and  is 
is  represented  by  the  equation 


TMtV/o) 


2  (40) 


C.   OBJECTIVES  OF  THE  PARAMETRIC  STUDY 

The  major  objective  of  the  present  analysis  was  to  refine  the  com- 
puter procedure  initiated  by  Knudson  so  that  a  more  comprehensive  para- 
metric optimization  could  be  carried  out.   The  primary  performance  in- 
dicators were  taken  to  be  the  specific  thrust  and  propulsive  efficiency. 
The  specific  thrust  and  propulsive  efficiency  of  the  water-augmented 
turbofan  were  compared  to  their  counterparts  for  a  "dry"  turbofan  which 
provided  a  thrust  ratio  and  an  efficiency  ratio  for  the  wet -to-dry  case. 
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The  system  parameters  described  below  were  varied  over  ranges  of  prac- 
tical values  to  determine  effects  on  thrust  and  efficiency  ratio. 

1 .  Fan  Total  Pressure  Ratio  and  Bypass  Ratio 

Perhaps  the  two  most  important  parameters  in  the  design  of  a 
turbofan  engine  are  the  fan  total  pressure  ratio  (FPR)  and  the  bypass 
ratio  (BR).   Consequently  it  was  expected  that  these  parameters  would 
also  be  of  importance  in  the  design  criteria  for  a  water-augmented 
turbofan  engine  and  that  information  concerning  the  effect  of  bypass 
ratio  and  fan  total  pressure  ratio  related  to  various  water-to-gas 
ratios  and  various  water  injection  velocities  could  be  very  useful  in 
optimizing  the  design  of  the  fan  section. 

2.  Water-to-Gas  Ratio 

The  effect  of  water-to-gas  ratio  variation  was  expected  to  be 
important  in  optimizing  the  thrust  and  efficiency  ratios  since  the 
basis  of  the  entire  concept  is  water-augmentation.   The  range  of  values 
of  water  volume  fraction  for  mist  flow  in  the  fan  duct  was  assumed  from 
zero  to  about  15  percent  at  standard  conditions  which  corresponds  to  a 
water-to-gas  ratio  range  of  zero  to  500. 

3 .  Mater  Injection  Velocity 

For  a  given  water-to-gas  ratio  a  wide  range  of  possible  water 

injection  velocities  (VTXTT)  exists.   The  effect  of  losses  caused  by  the 
J  INJ 

mixing  process  in  the  fan  mixing  duct  is  directly  related  to  the  veloc- 
ity ratio  between  the  water  and  air  in  the  injection  plane  making  the 
water  injection  velocity  a  significant  parameter.   The  forward  velocity 
of  a  vehicle  through  the  water  provides  only  a  certain  maximum  injection 
velocity  without  providing  additional  pump  work.   If  water  injection 
velocities  are  desired  which  are  greater  than  the  maximum  ram  velocity 
available,  allowing  for  inlet  duct  losses,  pump  work  must  be  added 
which  necessarily  reduces  the  efficiency  and  thrust  of  the  gas  genera- 
tor.  Only  if  the  advantages  of  increased  injection  velocity  outweigh  the 
losses  associated  with  the  pump  work  required  will  the  operation  be 
worthwhile.   Although  the  previous  work  by  Knudson  allowed  variation 
of  water  injection  velocity  the  pump  work  required  was  not  considered. 
One  of  the  objectives  of  the  present  analysis  was  to  provide  the  capabil- 
ity for  handling  injection  pump  work  as  required  by  variations  in  water 
injection  velocity. 
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4.  Optimum  Combinations  of  WGR  and  VT-_, 
_ IN  J 

The  efficient  operation  of  a  -water-augmented  turbofan  propul- 
sion system  would  require  a  knowledge  of  the  optumum  combinations  of 
water-to-gas  ratio  and  water  injection  velocity  to  provide  maximum 
thrust  at  a  given  ship  speed.   With  this  in  mind  it  was  decided  to 
provide  information  relating  the  optimum  combinations  required. 

5 .  Mixing  Duct  and  Nozzle  Thermal  and  Dynamic  Equilibrium 

The  water  and  air  enter  the  mixing  duct  at  different  velocities 
and  different  temperatures,  in  general.   During  the  mixing  process  the 
temperatures  and  velocities,  respectively,  tend  to  equalize  so  that 
different  ratios  of  water  temperature  to  air  temperature  and  water 
velocity  to  air  velocity  can  occur  at  the  mixing  duct  exit.   Similarly 
in  the  fan  nozzle  the  accelerations  and  temperature  changes  undergone 
by  the  two  phases  can  result  in  various  temperature  ratios  and  velocity 
ratios  at  the  nozzle  exit.   It  is  of  interest  to  determine  how  these 
ratios  affect  the  thrust  and  efficiency  of  the  propulsion  system. 

6 .  Fan  Duct  Pressure  Drop 

The  change  in  static  pressure  of  a  flow  in  a  duct  is  directly 
related  to  the  wall  shear  force  through  the  momentum  equation.   It  is 
possible,  therefore,  to  gain  some  understanding  of  the  effects  of  wall 
shear  friction  losses  in  two-phase  flow  by  looking  at  the  effects 
caused  by  variation  of  the  ratio  of  the  static  pressure  at  the  fan 
mixing  duct  exit  to  the  static  pressure  at  the  entrance  to  the  fan 
mixing  duct. 

D.   COMPUTER  PROGRAM  AND  MODIFICATIONS 

The  computer  program  used  in  the  present  study  is  based  upon  a 
standard  turbofan  analysis  which  was  modified  by  Knudson  to  consider 
the  injection  of  water  into  the  fan  discharge  duct.   An  explanation 
of  various  input  and  output  parameters  is  given  in  Appendix  I.   A 
copy  of  the  program  itself  and  a  sample  of  the  computer  output  are 
presented  in  Appendix  II. 

Originally  the  computer  program  consisted  of  two  separate  parts: 
a  "wet"  turbofan  program  and  a  "dry"  turbofan.   A  considerable  amount 
of  program  revision  was  required  to  enable  the  combination  of  the  two 
programs  so  that  the  ratios  of  water-augmented  thrust  and  efficiency 
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to  dry  thrust  and  efficiency,  respectively,  could  be  calculated  direct- 
ly.  Of  less  importance,  the  output  format  was  modified  to  present  the 
important  data  in  more  compact  form  to  allow  wider  ranges  of  the  var- 
ious parameters  to  be  covered  without  excessive  use  of  paper. 

1 .  Total  Entropy  Determination  at  Burner  Outlet 

In  the  original  program  the  calculation  used  in  determining 
the  total  entropy  at  the  burner  outlet  was  based  on  the  change  in  en- 
tropy from  the  burner  entrance  to  the  burner  exit  added  to  the  refer- 
ence entropy  of  the  atmosphere.   However,  the  use  of  the  reference 
entropy  of  the  atmosphere  in  Knudson's  analysis  was  incorrect  and  the 
reference  entropy  was  corrected  to  be  the  total  entropy  at  the 
burner  entrance . 

2 .  Temperature  Iterations 

Several  items  in  the  iterations  to  determine  the  mixing  duct 
air  temperature  and  the  fan  nozzle  air  temperature  were  improved  in 
accuracy  or  were  added  because  of  omissions  in  Knudson's  program. 
Of  these,  the  most  important  was  the  accuracy  indicator  used  to  deter- 
mine the  convergence  of  the  Newton-Raphson  iteration  for  air  tempera- 
ture at  both  the  duct  exit  and  the  nozzle  exit.   Knudson  examined  only 
the  numerator  of  the  convergence  term  whereas  the  proper  term  to  ex- 
amine is  the  absolute  value  of  the  entire  convergence  term.   This  was 
corrected  for  the  present  analysis  and  repeatable  results  were  ob- 
tained which  in  some  cases  differed  significantly  from  those  obtained 
by  Knudson. 

E.   RESULTS  AND  DISCUSSION 

The  significant  results  of  the  computer  study  are  presented  in 
graphic  form  in  Figs.  2  through  22.   For  all  results  shown  the  craft 
velocity  was  held  at  50  knots  although  the  computer  program  can  handle 
any  speed  range.   The  constant  velocity  allowed  more  meaningful  com- 
parisons to  be  made  among  the  results. 

The  following  parameters  were  also  held  constant  throughout  the 
analysis,  the  values  used  being  typical  of  presently  available  turbo- 
fan  engines.   The  compressor  total  pressure  ratio  was  maintained  at 
13.8  and  the  turbine  inlet  total  temperature  was  fixed  at  2500  degrees 
Rankine .   The  dry  fan  duct  total  pressure  ratio  was  taken  as  unity 
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since  no  separate  fan  duct  exists  on  a  dry  turbofan  engine.   Except 
for  the  study  of  duct  pressure  ratio  effects  the  wet  fan  duct  static 
pressure  ratio  was  arbitrarily  fixed  at  0.95  since  the  expected  two- 
phase  pressure  drop  was  not  known. 

The  maximum  available  velocity  for  water  injection  without  injec- 
tion pump  work  depends  upon  duct  and  piping  losses  in  transferring  the 
water  from  the  inlet  scoop  to  the  injection  nozzles.   The  velocity 
recovery  factor  chosen  for  the  present  analysis  was  0.8.   The  mach 
number  at  the  fan  outlet  (water  injection  plane)  was  assumed  to  be 
0.2   The  efficiencies  assigned  to  the  various  components  are  listed 
below  and  are  the  same  as  those  used  in  Refs.  2  and  4. 

Diffuser  0.95 

Compressor  0.90 

Burner  1.00 

Turbine  0.85 

Nozzle  0.98 

Fan  Diffuser  0.95 

Fan  0.90 

Fan  Nozzle  0.98 

Water  Injection  Pump      0.90 

1 .   Effects  of  Fan  Pressure  Ratio 

It  is  of  interest  to  determine  the  design  criteria  for  maximiz- 
ing thrust  ratio  with  respect  to  fan  pressure  ratio  and  bypass  ratio. 
The  variation  of  thrust  ratio  with  fan  pressure  ratio,  while  holding 
bypass  ratio  constant  at  4.0,  is  presented  for  various  water  injection 
velocities  in  Fig.  2  and  for  various  water-to-gas  ratios  in  Fig.  3. 
The  design  envelope  shown  in  both  graphs  is  determined  by  the  amount  of 
turbine  work  required  to  drive  the  compressor,  fan  and  water  injection 
pump.   The  turbine  work  available  is  dependent  on  the  maximum  pressure 
drop  allowable  across  the  turbine.   As  the  turbine  work  required  ex- 
ceeds the  turbine  work  available  the  turbine  exit  pressure  must  go  be- 
low atmospheric  pressure.   Since  this  is  impossible,  the  design  limit 
is  reached  when  the  turbine  exit  pressure  equals  atmospheric  pressure. 

It  is  observed  from  Fig.  2  that  for  a  given  water  injection 
velocity  or  water-to-gas  ratio  the  highest  possible  fan  pressure  ratio 
provides  the  maximum  thrust  ratio,  except  for  an  injection  velocity  of 
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100  knots,  below  a  fan  pressure  ratio  of  1.1.   The  increase  in  thrust 
ratio  for  low  values  of  fan  pressure  ratio  at  high  water  injection 
velocities  and  high  water-to-gas  ratios  is  considered  to  be  unreal- 
istic because  of  the  magnitudes  of  certain  of  the  parameters.   Specif- 
ically, a  fan  pressure  ratio  of  1.0  corresponds  to  a  duct  in  the  free 
stream  into  which  water  is  injected  with  no  losses  occurring  in  the 
fan.   Under  this  condition  higher  water  injection  velocities  reduce  the 
velocity  differential  between  the  water  and  the  faster  air  thus  poss- 
ibly reducing  mixing  losses  sufficiently  to  allow  an  increase  in 
thrust  ratio. 

The  optimum  fan  pressure  ratio  for  duct  and  nozzle  thermal 
and  dynamic  equilibrium,  respectively,  appears  to  be  in  the  range  1.4 
to  1.5  for  a  design  speed  of  50  knots.   This  range  gives  the  maximum 
thrust  ratio  while  still  allowing  the  widest  range  for  both  water  in- 
jection velocity  and  water-to-gas  ratio.   A  fan  pressure  ratio  of  1.5 
was  used  for  the  remainder  of  the  present  analysis.   The  effect  of 
variation  of  fan  pressure  ratio  on  efficiency  ratio  was  also  studied 
and  exhibited  trends  generally  similar  to  the  corresponding  thrust 
ratio  curves. 

2 .   Effects  of  Bypass  Ratio 

The  selection  of  the  optimum  bypass  ratio  is  also  an  important 
design  criteria  for  the  water-augmented  turbofan  engine.   The  variation 
of  thrust  ratio  with  bypass  ratio,  while  holding  fan  pressure  ratio 
constant  at  1.5,  is  presented  in  Fig.  4  for  various  water  injection 
velocities  and  Fig.  5  for  various  water-to-gas  ratios.   Again  the  de- 
sign envelope  reflects  the  maximum  available  turbine  work  limited  by 
turbine  exit  pressure  and  atmospheric  pressure  . 

Increases  in  both  injection  velocity  and  water-to-gas  ratio 
provide  increased  thrust  ratio  within  the  design  envelope.   Increases 
in  bypass  ratio  generally  increase  thrust  ratio  for  a  given  injection 
velocity  or  water-to-gas  ratio  except  for  injection  velocities  above 
about  70  knots  and  water-to-gas  ratios  greater  than  400.   The  decrease 
in  thrust  ratio  with  increasing  bypass  ratio  for  a  given  injection 
velocity  above  about  70  knots  (see  Fig.  4)  may  be  attributed  to  the 
fact  that  the  increase  in  pump  work  required  to  maintain  the  water- 
to  gas  ratio  constant  as  bypass  ratio  increases  detracts  more  from 
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the  thrust  of  the  gas  generator  than  the  amount  of  thrust  increase  of 
the  fan  section.   A  similar  explanation  applies  to  the  results  of 
Fig.  5  regarding  various  water-to-gas  ratios. 

The  optimum  design  bypass  ratio  for  maximizing  thrust  ratio 
at  a  design  speed  of  50  knots  appears  to  occur  in  the  range  of  bypass 
ratios  from  3.0  to  4.0.   This  range  allows  a  wide  selection  of  injec- 
tion velocities  and  water-to-gas  ratios  while  still  remaining  within 
the  design  envelope.   The  bypass  ratio  chosen  for  the  remainder  of  the 
present  analysis  was  4.0.   The  effect  of  variation  of  bypass  ratio  on 
efficiency  ratio  followed  trends  generally  similar  to  the  correspond- 
ing thrust  ratio  curves. 

3 .   Effects  of  Water-to-Gas  Ratio 

The  variation  of  thrust  ratio  with  water-to-gas  ratio  is  pre- 
sented in  Fig.  6  for  various  values  of  water  injection  velocity. 
Figure  7  shows  similar  information  with  respect  to  efficiency  ratio. 
It  is  apparent  that  both  water-to-gas  ratio  and  water  injection  veloc- 
ity greatly  affect  the  thrust  ratio  and  the  efficiency  ratio.   A  maxi- 
mum thrust  ratio  of  slightly  less  than  4.0  is  achieved  for  a  water-to- 
gas ratio  in  the  range  375  to  400  at  a  water  injection  velocity  slight- 
ly higher  than  60  knots  and  a  craft  velocity  of  50  knots  (see  Fig.  6). 
However  the  efficiency  ratio  shown  in  Fig.  7  continues  to  increase 
for  values  of  water-to-gas  ratio  exceeding  500.   At  a  water-to-gas 
ratio  of  500  the  maximum  efficiency  ratio  obtainable  is  about  6.2  for 
a  water  injection  velocity  slightly  below  60  knots  and  a  craft  velocity 
of  50  knots.   These  results  show  that  it  is  indeed  advantageous  to  in- 
crease the  injection  velocity  even  at  the  expense  of  the  pump  work 
required . 

The  operating  envelope  shown  in  Figs.  6  and  7  indicates  the 
maximum  turbine  work  available  for  driving  the  compressor,  fan  and 
water  injection  pump  and  arises  because  the  turbine  exit  pressure  is 
lower  than  atmospheric  pressure  at  the  higher  fan  pressure  ratios. 
The  explanation  given  earlier  in  Section  D.l  for  design  envelope  ap- 
plies to  the  operating  envelope  as  well. 

Each  water  injection  velocity  has  an  associated  water-to-gas 
ratio  at  which  thrust  ratio  is  maximized  and  another  water-to-gas  ratio 
for  which  efficiency  ratio  la  maximized.   This  may  be  explained  by  the 
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consideration  that  for  a  fixed  water  injection  velocity  the  initial 
velocity  ratio  between  the  water  droplets  and  the  air  is  constant.   As 
the  water-to-gas  ratio  increases  the  thrust  from  the  fan  nozzle  is  in- 
creased but  the  injection  pump  work  increase  required  causes  a  decrease 
in  the  thrust  of  the  gas  generator.   Initially  the  thrust  increase  from 
the  fan  section  predominates  but  eventually  the  thrust  decrease  from 
the  gas  generator  section  becomes  larger  than  the  thrust  gain  and  the 
net  thrust  begins  to  decrease. 
4.   Operating  Curves 

An  important  question  concerns  the  optimization  of  water  injec- 
tion velocity  and  water-to-gas  ratio  to  provide  maximum  thrust  ratio  at 
a  given  craft  velocity.  A  suggested  form  for  the  presentation  of  such 
information  is  shown  in  Fig.  8.  For  a  desired  craft  velocity  and  water 
injection  velocity  the  water-to-gas  ratio  required  to  provide  the  maxi- 
mum thrust  ratio  can  be  determined  from  the  curve. 

The  curves  do  not,  however,  provide  information  concerning  the 
optimum  combination  of  both  injection  velocity  and  water-to-gas  ratio 
required  to  provide  the  maximum  possible  thrust  ratio  at  a  given  craft 
velocity.   For  example,  at  a  water-to-gas  ratio  of  200  and  craft  veloc- 
ity of  50  knots  there  exist  two  injection  velocities,  41  knots  and  79 
knots,  respectively,  which  will  provide  a  maximum  thrust  ratio,  but  the 
greater  of  the  two  thrust  ratios  is  not  specified.  However,  comparison 
with  Fig.  6  indicates  thft  the  higher  thrust  ratio  coincides  with  the 
higher  injection  velocity.   Figure  7  also  shows  this  case  to  be  the  most 
efficient.   This  suggests  that  in  cases  where  two  possibilities  exist 
for  injection  velocity  the  higher  injection  velocity  should  be  used  to 
provide  the  higher  thrust  and  efficiency  ratio. 

The  importance  of  the  peaks  at  high  water-to-gas  ratios  on  the 
curves  for  craft  velocities  of  25  and  50  knots,  respectively,  lies  in 
the  fact  that,  for  water-to-gas  ratios  above  about  500,  the  fan  duct 
contains  such  a  high  percentage  of  water  that  the  water-jet  principle 
is  approached  with  its  increased  pumping  requirements  and  consequently 
increased  weight.   At  some  point  the  weight  increase  will  more  than 
offset  the  performance  advantage  and  the  water-augmented  turbofan  loses 
its  appeal.   For  this  reason  it  is  suggested  that  the  water-to-gas  ratios 
be  kept  as  low  as  possible. 
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5 .   Fan  Duct  and  Nozz le  Equilibrium  Parameters 

One  of  the  more  important  aspects  of  the  present  computer  anal- 
ysis has  been  the  determination  of  the  effect  of  fan  duct  and  nozzle 
equilibrium  parameters  upon  the  thrust  ratio  and  the  efficiency  ratio. 
The  important  equilibrium  parameters  are  taken  to  be  the  fan  duct  exit 

velocity  ratio,  V\,„/VA„,  the  fan  duct  exit  temperature  ratio,  T  „/T  „, 
'   W8   A8  '   W8   A8 

fan  nozzle  exit  velocity  ratio,  Vr7r./V._,  and  the  fan  nozzle  exit  tem- 

J  '   W9   A9' 

perature  ratio,  T  Q/T   ;  these  ratios  are  referred  to  in  the  succeeding 
wy  Ay 

test  as,  respectively,  duct  velocity  ratio,  duct  temperature  ratio, 
nozzle  velocity  ratio  and  nozzle  temperature  ratio  and  are  water-to- 
air  ratios. 

These  parameters  are  important  for  several  reasons.   First, 
little  is  known  about  the  nature  of  the  mixing  process;  consequently, 
analysis  must  be  made  on  an  empirical  basis.   Second,  Lockhart  and 
Martinelli  [Ref.  61  and  Chenoweth  and  Martin  |Ref.  7J  indicate  that 
mixing  duct  pressure  losses  are  greater  for  the  flow  of  a  two-phase 
mixture  than  for  either  phase  flowing  separately.   The  extent  of  the 
increase  and  the  factors  which  govern  the  increase  are  not  very  well 
known.   Further  discussion  of  this  subject  is  presented  in  Section  II 
of  this  paper. 

a.  Fan  Duct  Velocity  Ratio 

The  variation  of  thrust  ratio  and  efficiency  ratio  with 
duct  velocity  ratio  is  presented  in  Figs.  9  and  10,  respectively,  for 
various  values  of  water  injection  velocity.   Under  the  restriction  that 
all  other  equilibrium  parmeters  have  the  value  unity,  the  duct  velocity 
ratio  has  only  a  very  slight  effect  on  thrust  ratio  and  efficiency  ratio 
The  duct  velocity  ratio,  therefore,  does  not  appear  to  be  a  significant 
parameter  affecting  thrust  ratio  and  efficiency  ratio. 

b.  Fan  Duct  Temperature  Ratio 

The  effects  of  duct  temperature  ratio  on  thrust  ratio  and 
efficiency  ratio  are  shown  in  Figs.  11  and  12,  respectively,  for  various 
water  injection  velocities.   The  initial  conclusion  from  Fig.  11  is  that 
it  would  be  highly  desirable  to  maintain  the  duct  temperature  ratio  as 
low  as  possible  to  attain  the  maximum  thrust  ratio.   However,  further 
consideration  of  the  restrictions  imposed,  namely  duct  velocity  equi- 
librium, nozzle  velocity  equilibrium  and  nozzle  temperature  equilibrium, 
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create  unrealistic,  if  not  impossible,  conditions  for  the  flow  to  meet. 

First,  it  is  extremely  unlikely  that  a  low  value  of  duct 
temperature  ratio  could,  at  the  same  time,  produce  equilibrium  of  the 
other  three  parameters.   Second,  the  experimental  results  of  Reese  and 
Richard  JRef .  8]  indicate  that  thermal  equilibrium  occurs  in  two-phase 
flows  with  contact  times  as  low  as  the  order  of  one  millisecond.   A  typ- 
ical expected  contact  time  of  the  order  of  20  milliseconds  suggests 
thermal  equilibrium  at  all  times  in  the  fan  duct  exit. 

A  further  consideration  is  that  for  low  values  of  duct  tem- 
perature ratio  the  bulk  of  the  mixing  process  must  occur  in  the  nozzle. 
Since  the  calculations  through  the  nozzle  include  the  assumption  of  no 
net  entropy  change,  the  thrust  ratio  resulting  is  larger  than  that  which 
could  realistically-  occur.   These  reasons  suggest  that  the  results  given 
in  Figs.  11  and  12  do  not  give  an  accurate  indication  of  realistic  ef- 
fects . 

c.   Fan  Nozzle  Velocity  Ratio 

Figures  13  and  14,  respectively,  present  the  effect  on 
thrust  ratio  and  efficiency  ratio  of  varying  the  ratio  of  water  velocity 
to  air  velocity  at  the  nozzle  exit  for  various  water  injection  veloci- 
ties.  For  a  given  injection  velocity  neither  thrust  ratio  nor  efficiency 
ratio  varies  a  great  amount,  however,  the  variation  is  greater  than  the 
similar  case  for  variation  of  duct  velocity  ratio.   For  an  injection 
velocity  of  100  knots  the  thrust  ratio  increases  from  about  3.13  at  a 
velocity  ratio  of  0.5  to  about  3.22  at  a  velocity  ratio  of  1.0  and  the 
efficiency  ratio  increases  from  about  3.77  to  about  3.88.   Both  curves 
become  nearly  horizontal  near  a  duct  velocity  ratio  of  1.0. 

It  should  be  pointed  out  that  the  restriction  of  equilib- 
rium conditions  for  the  other  parameters  is  not  so  unrealistic  at  the 
nozzle  exit  as  it  was  at  the  duct  exit.   This  is  because  the  flow  is 
accelerated  through  the  nozzle  and  the  air,  because  of  its  lower  density, 
accelerates  more  quickly  than  the  water  causing  a  velocity  difference  at 
the  nozzle  exit  regardless  of  the  equilibrium  state  at  the  fan  duct  exit 
(nozzle  entrance).   As  was  the  case  with  duct  velocity  ratio,  the  nozzle 
velocity  ratio  appears  to  have  little  significant  effect  upon  thrust 
ratio  and  efficiency  ratio. 
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d.   Fan  Nozzle  Temperature  Ratio 

The  variation  of  thrust  ratio  versus  nozzle  temperature 
ratio  is  presented  for  various  water  injection  velocities  in  Fig.  15, 
and  for  various  water-to-gas  ratios  in  Fig.  16.   In  both  cases  the 
thrust  ratio  decreases  smoothly  with  increasing  nozzle  temperature  ratio 
indicating  that  the  most  desirable  operating  range  is  near  thermal 
equilibrium  where  the  curves  are  nearly  flat.   Again,  based  on  the  re- 
sults of  Reese  and  Richard,  it  is  probable  that  the  nozzle  temperature 
ratio  will  be  near  unity  in  spite  of  the  temperature  changes  of  the  air 
through  the  nozzle. 

The  variation  of  efficiency  ratio  with  nozzle  temperature 
ratio  is  shown  in  Figs.  17  and  18,  respectively,  for  various  water  in- 
jection velocities  and  various  water-to-gas  ratios.   The  results  tend 
to  follow  trends  similar  to  those  for  the  thrust  ratio  except  for  the 
higher  values  of  water  injection  velocity  and  water-to-gas  ratio.   For 
a  water  injection  velocity  of  about  80  knots  in  Fig.  17  the  efficiency 
ratio  is  relatively  unchanged  over  the  range  of  nozzle  temperature 
ratios.   An  injection  velocity  of  100  knots  shows  increasing  efficiency 
ratio  with  increasing  nozzle  temperature  ratio.   Similarly  a  water-to- 
gas ratio  of  400  or  higher  in  Fig.  18  exhibits  a  similar  trend  of  in- 
creasing efficiency  ratio. 

The  information  for  temperature  ratios  significantly  greater 
than  unity  is  considered  to  be  unrealistic  in  view  of  the  predicted  short 
contact  times  required  to  achieve  thermal  equilibrium.   In  addition  the 
large  pressure  drop  through  the  nozzle,  which  would  be  required  for  such 
a  great  air  temperature  change,  is  highly  unlikely  to  occur. 
6 .   Effects  of  Fan  Duct  Pressure  Ratio 

An  important  parmeter  affecting  the  performance  of  the  water- 
augmented  turbofan  engine  is  the  pressure  loss  in  the  fan  mixing  duct. 
Since  the  wall  shear  force  in  two-phase  flow  is  an  unknown  quantity 
which  has  not  been  treated  analytically  in  this  paper,  a  range  of  values 
for  the  ratio  of  static  pressure  at  the  mixing  duct  exit  to  the  static 
pressure  at  the  duct  entrance  has  been  used  to  determine  the  effect  on 
thrust  ratio  and  efficiency  ratio.   Figures  19  and  20  present  thrust 
ratio  versus  duct  pressure  ratio  for  various  water  injection  velocities 
and  various  water-to-gas  ratios,  respectively.   Figures  21  and  22  pre- 
sent parallel  results  for  efficiency  ratio,  respectively. 
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It  is  theoretically  possible  for  the  static  pressure  to  in- 
crease in  the  direction  of  flow  if  large  velocity  variations  exist  over 
the  injection  plane  and  the  wall  shear  stress  is  small.   Thus,  these 
plots  have  been  extended  to  pressure  ratios  greater  than  one.   It  is 
easily  seen  that  the  most  desirable  range  of  operation  to  provide  max- 
imum thrust  ratio  is  at  the  maximum  duct  pressure  ratio.   This  is 
simply  the  expected  result  that  minimum  duct  friction  provides  the 
best  operating  conditions.   Also,  as  predicted  by  a  momentum  analysis, 
the  highest  water  injection  velocity  provides  the  least  sensitivity  to 
duct  pressure  losses  as  well  as  providing  the  best  thrust  ratio. 

The  effect  of  duct  pressure  ratio  on  efficiency  ratio  is  simi- 
lar to  the  effect  on  thrust  ratio;  however,  for  a  fixed  water-to-gas 
ratio  the  efficiency  ratio  is  virtually  unchanged  with  varying  duct 
pressure  ratio.   Again  increasing  water  injection  velocity  decreases 
the  effect  on  efficiency  ratio  of  varying  duct  pressure  ratio.   As  the 
water  injection  velocity  increases  to  about  80  knots,  the  injection 
pump  work  required  begins  to  predominate  over  the  efficiency  gain  af- 
forded by  the  high  water  injection  velocity,  causing  the  efficiency 
ratio  to  decrease  with  increasing  duct  pressure  ratio  (see  Fig.  21). 

F.   CONCLUSIONS 

The  injection  of  water  into  the  fan  discharge  duct  of  a  turbofan 
engine  provides  a  feasible  propulsion  system  for  high  speed  sea-borne 
vessels.   Among  the  most  important  parameters  affecting  the  thrust 
and  propulsive  efficiency  of  water-augmented  turbofans  are  the  water- 
to-gas  ratio  and  the  water  injection  velocity.   Each  water  injection 
velocity  has  an  associated  water-to-gas  ratio  which  maximizes  thrust 
and  another  water-to-gas  ratio  which  maximizes  propulsive  efficiency. 
Thrust  ratios  and  efficiency  ratios  of  the  water-augmented  turbofan 
to  the  dry  turbofan  of  three  or  four  are  easily  attainable.   It  should 
be  noted  that  these  figures  do  not  account  for  the  increased  drag  that 
would  be  caused  by  the  water  inlet  scoop. 

Velocity  slip  ratio  in  either  the  fan  mixing  duct  or  the  fan  nozzle 
has  negligible  effect  upon  the  thrust  ratio  and  efficiency  ratio.   The 
effect  of  temperature  difference  between  water  and  air  is  somewhat 
greater  than  for  the  velocity  ratio.   However,  the  likelihood  of  non- 
equilibrium  for  either  fan  duct  or  fan  nozzle  temperature  ratio  is 


31 


slight  due  to  the  short  contact  times  required  for  thermal  equilibrium. 

Thrust  ratio  and  efficiency  ratio  both  increase  with  increasing 
fan  duct  pressure  ratio.   Fan  duct  pressure  ratio  can  be  related  to 
fan  duct  wall  friction  forces  through  the  momentum  equation,  indicating 
that  minimum  wall  friction  is  desirable.   Thrust  ratio  and  efficiency 
ratio  become  less  sensitive  to  fan  duct  pressure  ratio  with  increasing 
water  injection  velocity  and  increasing  water-to-gas  ratio.   Since 
little  is  known  about  two-phase  pressure  drop,  experimental  work  should 
be  done  to  increase  the  knowledge  in  the  field.   Section  II  of  this 
paper  describes  the  development  of  a  test  facility  for  wall  shear  force 
and  pressure  drop  measurements  in  two-phase  flow. 

For  a  design  speed  of  50  knots,  the  ideal  fan  total  pressure  ratio 
falls  within  the  range  1.4  to  1.5,  if  use  of  a  wide  range  of  water 
injection  velocities  and  water-to-gas  ratios  is  desired.   At  the  same 
design  speed,  the  optimum  fan  bypass  ratio  occurs  in  the  range  3.0  to 
4.0,  to  allow  use   of  a  wide  range  of  water  injection  velocities  and 
water-to-gas  ratios. 

For  any  given  water  injection  velocity  a  set  of  curves  may  be  con- 
structed which  provide  the  water-to-gas  ratio  required  to  achieve  maxi- 
mum thrust  for  a  specified  craft  velocity.   The  present  analysis  has 
not  considered  the  relationship  between  increasing  the  water  injection 
velocity  and  the  resulting  weight  increase  due  to  pump  requirements. 
As  the  water  injection  velocity  and  water-to-gas  ratio  increase  the 
pure  water-jet  concept  is  approached  and  the  performance  advantages 
of  the  water-augmented  turbofan  concept,  primarily  the  low  system 
weight,  are  lost.   A  further  study  is  thus  required  to  determine  the 
point  at  which  weight  and  cost  increases  for  water  injection  pumps 
begin  to  detract  significantly  from  the  overall  system  effectiveness. 
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II.   DEVELOPMENT  OF  A  TEST  FACILITY  FOR  TWO- PHASE  FLOW 

A.   TWO- PHASE  PRESSURE  DROP  PREDICTION 

The  ability  to  predict  the  pressure  drop  occurring  in  a  two-phase 
duct  flow  is  vital  to  the  accurate  analysis  of  the  water-augmented 
turbofan  concept  discussed  in  Section  I.   To  the  author's  knowledge  no 
complete,  accurate,  theoretical  analysis  of  two-phase  pressure  drop 
has  as  yet  been  proposed.   However,  Dukler,  Wicks  and  Cleveland  IRef.  91 
and  other  available  literature  indicate  that  two  empirical  data  correl- 
ations for  the  prediction  of  two-phase  pressure  drop  have  given  con- 
sistently better  results  than  other  correlations.   These  correlations 
were  developed,  respectively,  by  Lockhart  and  Martinelli  [Ref.  6j   and 
Chenoweth  and  Martin  [Ref .  7J   .   An  explanation  of  each  of  the  correl- 
ations and  some  of  the  special  assumptions  made  in  them  is  presented 
in  the  following  text. 

1 .   Lockhart-Martinelli  Two-Phase  Flow  Correlation 

The  correlation  of  Lockhart  and  Martinelli  is  dependent  upon 
which  of  four  types  of  flow  exists  during  the  simultaneous  flow  of  a 
liquid  and  a  gas  (or  vapor).   These  flow  regimes  are: 

1.  Flow  of  both  the  liquid  and  the  gas  may  be  turbulent 

(turbulent- turbulent  flow) 

2.  Flow  of  the  liquid  may  be  viscous  and  flow  of  the  gas 
may  be  turbulent  (viscous- turbulent  flow) 

3.  Flow  of  the  liquid  may  be  turbulent  and  flow  of  the  gas 
may  be  viscous  (turbulent-viscous  flow) 

4.  Flow  of  both  the  liquid  and  the  gas  may  be  viscous 
(viscous-viscous  flow) 

"Viscous"  flow  is  the  term  used  by  Lockhart  and  Martinelli  for  what 

is  more  properly  called  laminar  flow. 

The  assumptions  made  are  that  the  "static  pressure  drop  for 

the  liquid  phase  must  equal  the  static  pressure  drop  for  the  gaseous 

phase  regardless  of  the  flow  pattern,  as  long  as  an  appreciable  radial 

static  pressure  difference  does  not  exist"  and  "the  volume  occupied 

by  the  liquid  plus  the  volume  occupied  by  the  gas  at  any  instant  must 

equal  the  total  volume  of  the  pipe."   It  is  further  assumed  that  the 

transition  from  laminar  to  turbulent  flow  occurs  in  the  range  of 

Reynolds  numbers  between  1000  and  2000.   The  Reynolds  number  is  calcu- 
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lated  using  the  mass  flow  rate  of  the  phase  in  question  and  the  inside 
diameter  of  the  pipe.   For  example,  the  liquid  phase  Reynolds  number  is 
determined  from  the  equation 


ReL  = 


4-m, 


(41) 


*  OyUL 

where  m,  is  the  liquid  mass  flow  rate,  in  pounds-mass  per  second, 
D  is  the  inside  pipe  diameter  in  feet  and  yty   is  the  absolute  viscosity 
of  the  liquid  in  pounds-mass  per  foot-second. 

Having  determined  which  of  the  flow  regimes  is  appropriate 
from  the  Reynolds  number  calculation  for  each  phase,  the  all-liquid 
and  all-gas  pressure  drops,  respectively,  are  calculated.   The  all- 
liquid  pressure  drop  in  a  pipe  L  feet  in  length  is  based  upon  the 
liquid  flowing  at  a  rate  itl  with  a  density  P,     as  in  the  equation 

AK^fiK-'fri  (42) 

where  the  liquid  velocity,  V  ,  is  determined  from  continuity,  assuming 

J-j 

that  the  liquid  occupies  the  entire  pipe,  and  f   is  the  liquid  friction 

Jj 

factor  based  on  Re  .   Similarly,  an  all-gas  pressure  drop  is  found  from 
J-j 

APe    -.■/*&£■£  (43) 

The  ratio  of  the  all-liquid  pressure  drop  to  the  all-gas  pres- 
sure drop  is  the  square  of  the  two-phase  flow  modulus,  "JC  ,    or 


X.   ^ 


APl 


-J/o 


A  Pa 


(44) 


am- 


Figure  23  shows  the  relationship  between  ^C  and  the  correlation  par 
eter  Q    for  each  of  the  various  flow  regimes,  as  determined  experiment- 
ally by  Lockhart  and  Martinelli.   If  (g.     is  taken  from  the  curve  the 
two-phase  pressure  drop  is 

APrp    =  $?&Pl.  (45) 

Similarly,  if  $c  is  taken  from  the  curve  the  two-phase  pressure  drop  is 


APtp  s  &apg 


(46) 


It  should  be  noted  that  Lockhart  and  Martinelli  state  that 
more  data  are  needed  to  establish  the  validity  of  their  correlation  at 
very  high  and  very  low  values  of  the  flow  modulus  31   (nearly  all- 
liquid  and  nearly  all-gas  flow,  respectively). 
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2 .   Chenoweth-Martin  Two-Phase  Flow  Correlation 

Chenoweth  and  Martin  developed  an  "improved  correlation  for 
two-phase  pressure  drop  in  horizontal  pipes"'  which  is  especially  valid 
for  pressures  up  to  100  psia,   The  correlations  can  be  used  for  any 
two-phase  mixture  as  long  as  the  flow  is  turbulent  and  it.  predicts 
single-phase  values  when  the  flow  is  all  liquid  or  all  gas.   The  two 
input  parameters  for  the  use  of  the  correlation  are  the  liquid  volume 
fraction  (LVF)  and  the  ratio  of  a  fictitious  all-gas  pressure  drop  to 
a  fictitious  all-liquid  pressure  drop. 

The  liquid  volume  fraction  is  calculated  from  the  flow  rates 
and  densities  of  the  two-phases  using  the  equation 

Volume  flow  rate  of  liquid       U'Pu  * 


LVF      = 


Total  volume    flow  rate  >ntJ  f  "***&/&         1  4-   "H&MMit 


The  ratio  of  the  fictitious  pressure  drops  is  determined  by  computing 
the  value  of 

where 


•"kfa  (47) 
pressure   drops    is   determined   b) 

ft  Kfa  ft* 

1*    H  l*'t*2K  (48) 

and 

%*  *  <i*.kt.EK  (49) 

The  subscripts  L*  and  G*   denote  the  fictitious  all-liquid  and  all-gas 
states,  respectively.   The  K' s  are  loss  coefficients  for  valves  and 
fittings.   The  fictitious  friction  factor  is  determined  in  the  normal 
manner  using  an  artificial  Reynolds  number  based  on  the  total  flow  rate 
of  liquid  and  gas  and  the  physical  properties  of  the  phase  in  question. 
For  example,  f    is  determined  using  the  Reynolds  number 

Li  'v 


7T  D^(L 

The  fictitious  all-liquid  pressure  drop  is  then  computed  as 

where  VT .,  is  a  fictitious  velocity  determined  from  the  equation 
J_i 


(51) 


tf«    = 


Tf\t_  -f  7r>fi 


^     '  /Oj_A  (52) 

Figure    24  can   now  be   used    to   predict    the    two-phase    pressure    drop,  A  Prp , 
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In  contrast  to  the  Lockhart-Martinelli  correlation,  the 
Chenoweth-Martin  correlation  covers  only  all- turbulent  two-phase 
mixtures  but  ties  into  single-phase  values  at  either  end.   It  is 
also  valid  over  a  wider  range  of  pressures. 

B.   EQUIPMENT  DESIGN  OBJECTIVES  AND  DESCRIPTION 

A  test  rig  was  developed  for  the  experimental  analysis  of  two-phase 
flows,  in  particular  the  mixing  process  which  immediately  follows 
injection  of  water  into  the  air.   Initially,  fully  established  two- 
phase  flows  will  be  investigated  by  directly  measuring  the  wall  fric- 
tion force  and  the  overall  pressure  drop.   Subsequent  investigation 
will  concentrate  on  determination  of  the  pressure  drop  due  to  the  mixing 
process  in  a  developing  flow.   The  direct  measurement  of  the  wall  fric- 
tion force  will  allow  separation  of  the  pressure  drop  due  to  mixing  from 
the  overall  pressure  drop.   General  views  of  the  apparatus  are  shown 
in  Figures  25  and  26.   A  flow  diagram  for  the  initial  tests  is  shown 
in  Figure  27 . 

Air  enters  the  system  through  a  2-inch  pipe  from  a  plenum  fed  by  a 
reciprocating  compressor  capable  of  supplying  air  at  pressures  up  to 
150  psig  and  flow  rates  up  to  700  scfm.   Water  is  supplied  to  the  sys- 
tem through  a  3-inch  pipe  from  a  centrifugal  pump  capable  of  supplying 
up  to  500  gpm  at  80  psig.   The  mass  flow  rates  of  the  air  and  water  are 
accurately  determined  using  thin-plate  orifice  meters.   The  air  and 
water  flows  are  united  by  joining  the  two  pipes  into  a  3-inch  pipe. 
The  two-phase  mixture  flows  together  for  about  28  feet  becoming  fully 
established,  then  flows  through  a  "floating"  test  section  10  feet  long 
in  which  the  wall  shear  force  is  measured.   The  downstream  end  of  the 
test  section  exhausts  to  the  atmosphere  where  the  water  is  collected 
and  returned  to  the  water  pump. 

1 .   Flow  Orifice  Details 

The  flow  measurement  orifices  shown  in  Fig.  28  were  designed 
according  to  standard  ASME  specifications   [Ref.  1 0 1  .   Figures  29A 
and  29B  give  the  dimensions  of  the  orifices  and  the  placement  of  the 
pressure  taps  with  respect  to  the  pipe.   Each  orifice  plate  was 
machined  from  type  304  stainless  steel  plate  for  which  a  coefficient 
of  expansion  curve  is  available  in  Ref.  10.   Each  of  the  flanges  con- 
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taining  an  orifice  was  machined  with  a  circular  groove  centered  with 
respect  to  the  pipe  axis  to  insure  accurate  alignment  of  the  orifice 
with  the  pipe  axis. 

2 .  Initial  Mixing  of  Air  and  Water 

The  mixing  of  the  air  and  water  is  accomplished  by  merely 
joining  the  two  pipes  at  a  45-degree  angle  (see  Fig.  30).   The  mixing 
process  occurs  solely  due  to  the  turbulent  nature  of  the  flow;  if  the 
liquid  volume  fraction  is  sufficiently  small  a  mist-flow  will  result. 
It  is  the  mist-flow  regime  (liquid  volume  fractions  below  about  15 
percent)  which  is  directly  applicable  to  the  water-augmented  turbofan 
engine  described  in  Section  I.   However,  the  apparatus  is  designed  to 
handle  as  wide  a  range  of  water-to-air  mass  flow  ratios  as  permitted 
by  the  available  water  pump  and  air  compressor  capacities. 

3 .  Test  Section  Details 

The  test  section  consists  of  a  10-foot  long,  3.106-inch  steel 
pipe  which  is  supported  two  feet  from  each  end.   Each  support  has  four, 
one-eighth  inch  diameter  stainless  steel  rods,  one  vertically  above,  one 
vertically  below,  and  one  horizontally  on  each  side  as  shown  in  Fig.  31. 
The  rods  are  mounted  on  bearings  which  allow  the  test  section  to  "float" 
longitudinally . 

The  wall  shear  force  measurement  is  made  through  two  cantilever 
flexures  mounted  horizontally,  one  on  each  side  of  the  test  section 
(see  Fig.  32).   Each  flexure  has  four  strain  gages  mounted  to  detect 
changes  in  bending  moment  caused  by  the  force  exer'ed  by  the  test  sec- 
tion.  Figure  33  shows  the  wiring  diagram  for  each  flexure. 

The  labyrinth  seal  shown  in  Fig.  34  was  designed  to  provide  a 
pressure  seal  at  the  inlet  to  the  test  section  which  does  not  restrict 
longitudinal  motion.   The  dimensional  details  of  the  labyrinth  are 
given  in  Fig.  35.   Experimental  determination  of  the  leakage  rate  is 
discussed  in  Section  C  below. 

The  pressure  at  the  inlet  to  the  test  section  is  measured  to 
allow  the  determination  of  the  pressure  drop  through  the  test  section. 
Additionally  the  test  section  inlet  pressure  is  required  in  order  to 
calculate  the  force  on  the  labyrinth  face  at  the  test  section  inlet. 
This  pressure  force  must  be  subtracted  from  the  total  force  indicated 
on  the  flexures  to  obtain  the  true  wall  shear  force. 
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C   CALIBRATION  PROCEDURE 

It  is  necessary  to  determine  the  amount  of  air  leakage  which  occurs 
at  the  labyrinth  seal  over  the  range  of  expected  internal  to  external 
pressure  ratios.   This  was  to  have  been  accomplished  by  sealing  the 
outlet  of  the  test  section,  varying  the  internal  pressure  and  measuring 
the  mass  flow  rate  of  air  leaving  through  the  labyrinth  seal.   However, 
excessive  leakage  in  the  air  supply  line  prevented  the  accurate  meas- 
urement of  the  leakage  rate  through  the  labyrinth  seal  and  insufficient 
time  remained  to  complete  the  necessary  repairs. 

The  flexures  on  which  the  strain  gages  are  mounted  require  calibra- 
tion before  and  after  each  series  of  test  runs.   The  calibration  proced- 
ure for  the  flexures  involves  placing  a  known  force  on  the  test  section 
using  a  system  of  weights  and  a  pulley  as  shown  in  Fig.  36.   A  typical 
calibration  curve  of  force  versus  strain  gage  output  is  presented  in 
Fig.  37. 

D.   DATA  COLLECTION  AND  REDUCTION  PROCEDURE 

The  important  variables  for  the  operation  of  the  two-phase  flow 
test  rig  are  the  mass  flow  rates  of  air  and  water,  the  force  exerted 
on  the  test  section  due  to  wall  friction  and  labyrinth  pressure  effects, 
and  the  pressure  loss  in  the  test  section.   To  determine  the  mass  flow 
rates  it  is  necessary  to  record  the  static  pressure  upstream  of  the  air 
orifice,  the  pressure  drop  across  both  orifices  for  both  the  flange 
taps  and  the  D  and  half-D  taps,  respectively,  the  temperature  of  each 
fluid  and  the  atmospheric  temperature  at  each  manometer. 

The  calculation  of  the  two-phase  pressure  drop  predictions  of  both 
Lockhart  and  Martinelli  and  Chenoweth  and  Martin  require  that  the  inlet 
static  pressure  in  the  test  section  and  the  temperature  of  the  air-water 
mixture  at  the  outlet  of  the  test  section  be  recorded.   Atmospheric 
pressure  must  be  recorded  for  use  in  mass  flow  rate  calculations.   The 
temperature  is  recorded  at  the  outlet  of  the  test  section  because  the 
temperature  probe  will  disturb  the  flow  if  it  is  installed  at  the  inlet 
to  the  test  section.   It  should  be  mentioned  at  this  point  that  the 
temperature  measured  in  a  two-phase  mixture  can  be  subject  to  discussion 
because  it  is  not  known  what  effect  the  water  droplets  impinging  on  the 
probe  have  on  the  indicated  temperature.   The  accuracy  of  the  pressure 
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drop  predictions  depends  heavily  upon  an  accurate  temperature  of 
the  flow. 

A  computer  program  converts  the  raw  experimental  data  into  mass 
flow  rates  and  two-phase  pressure  drop  predictions.   The  input  and 
output  variables  used  for  the  computer  program  are  defined  in  Appendix 
III  and  the  computer  program  itself  along  with  a  sample  output  is  pre- 
sented in  Appendix  IV. 

The  pressure  drop  prediction  of  Chenoweth  and  Martin  is  not  com- 
pletely carried  out  by  the  computer  program  because  of  the  difficulty 
in  obtaining  an  analytical  expression  for  the  curves  of  the  correlation 
in  Fig.  24.   The  program  does  supply  the  variables  necessary  to  enter 
Fig.  24.   The  value  of  APf-p/Ap^x         thus  found  is  multiplied  by  the 
fictitious  all-liquid  pressure  drop,   AP/jt      ,    which  is  also  computed 
in  the  program,  to  obtain  the  predicted  two-phase  pressure  drop. 

E.   CURRENT  STATUS  AND  SUGGESTIONS  FOR  FURTHER  WORK 

The  work  reported  herein  concerning  the  development  of  the  two- 
phase  flow  test  rig  suffered  several  setbacks  due  to  difficulty  in 
obtaining  parts.   Because  of  these  delays  the  installation  of  the 
water  return  system  has  not  been  completed. 

Several  suggestions  for  further  work  have  become  apparent  during 
the  course  of  the  development  of  the  two-phase  flow  test  rig, 

1 .  Labyrinth  Air  Leakage  Determination 

It  will  be  necessary  to  determine  the  leakage  rate  of  air 
through  the  labyrinth  seal  prior  to  taking  experimental  data.   It 
should  be  noted  that  the  possibility  exists  that  the  present  metering 
system  may  not  be  adequate  for  measuring  minimum  flow  rates  of  both 
air  and  water. 

2 .  Labyrinth  Water  Leakage.  Determination 

During  operation  with  two-phase  or  all-liquid  flow  some  water 
will  probably  leak  out  the  labyrinth  seal  and  provision  must  be  made 
to  measure  this  leakage. 

3 .  Computerization  of  Friction  Factor  Curve 

The  friction  factors  used  in  calculating  pressure  drops  in  the 
computer  program  solutions  for  predicted  two-phase  pressure  drop  were 
assumed  to  be  constant  values,   This  is  true  only  in  the  region  of 
fully  turbulent  flow.   The  value  of  this  friction  factor  used  for  high 


39 


Reynolds  numbers  should  first  be  checked  by  simple  pressure  drop  tests. 
The  computer  program  must  then  be  improved  by  the  computerization  of 
the  friction  factor  curve  so  that  the  laminar  flow  and  transition  re- 
gions will  be  included. 

4 .  Computerization  of  Chenoweth-Martin  Correlation  Curve 

It  would  be  helpful  if  an  analytical  expression  could  be  found 
for  the  Chenoweth  and  Martin  correlation  curves  of  Fig.  24. 

5 .  Evaluation  of  Predicted  Pressure  Drop 

The  knowledge  gained  from  the  continuation  of  the  work  reported 
herein  should  be  used  to  evaluate  the  pressure  drop  correlations  of 
Lockhart  and  Martinelli  and  Chenoweth  and  Martin.   Also,  a  prediction 
for  the  two-phase  wall  shear  has  very  recently  been  proposed  by  Muench 
and  Ford  [Ref.  Ill  and  should  be  evaluated.   The  analysis  of  the  two- 
phase  flow  calculations  in  the  water-augmented  turbofan  should  be  re- 
examined in  the  light  of  experimental  results. 

6 .  Pressure  Drop  due  to  Two-Phase  Mixing 

The  present  test  rig  is  set  up  to  obtain  fully  developed  two- 
phase  flow  in  the  test  section.   However  if  the  water  is  injected  im- 
mediately upstream  from  the  test  section  it  should  be  possible  to 
separate  the  pressure  drop  into  that  caused  by  the  mixing  of  the  two 
phases  and  that  caused  by  wall  friction.   This  knowledge  would  be  use- 
ful in  further  analysis  of  injection  systems. 
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FIGURE 24.  CORRELATION   OF  CHENOWETH    AND  MARTIN    [REF  6]    FOR   TURBULENT 
TWO-PHASE    PRESSURE   DROP    IN   HORIZONTAL   PIPES 
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A  -Flange  Taps 
B  -  D  and  I/2D  Taps 


Air  Pipe 

D=  2.087 

D/2= 1.043" 
d  =1.460 


D  =   3.086 

D/2  =1.543" 
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FIGURE       29  A. 

GENERAL  ARRANGEMENT  OF    FLOW  ORIFICE 
INSTALLATIONS 


SECTION   X-X 


FIGURE      23  B. 
LOCATION  OF  ORIFICE     PRESSURE    TAPS  AROUND  PIPE 
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FIGURE  31 
TEST  SECTION  SUPPORTS 
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Siqnol  (-) 


A.  Switching  And    Boloncing  Unit 

B.  Daytronic    Strain  Gage  Digital  Indicator 

Model  700 

C.  Compression    Strain  Gages 

T.    Tension    Strain  Gages 

Strain  Gages  ••     Budd    120  £1    Foil  Gages 
Bonding  Agent:  Eastman  9K) 
Waterproofing:     Dow  Corning   RTV 

FIGURE    33 
WIRING    DIAGRAM      FOR   EACH   FLEXURE 
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FIGURE  36 
FLEXURE  CALIBRATION  SET-UP 
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FIGURE  37 
TYPICAL  FLEXURE  CALIBRATION    CURVE 
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APPENDIX  I 


TURBOFAN  PROGRAM  NOW     ,     fRE 


Input  Parameters 
Card  1: 


ETAD 

ETAC 

ETAB 

ETAT 

ETAN 

ETAFD 

ETAF 

ETAFN 

EPUMP 


Card  2: 


DSO 
PSO 

TSO 
HVF 
TREFA 


Card  3: 


AMN7 
TWO 
PTR32 
PTR87 


Card  4: 


NVO 

NCPR 

NFPR 

NTT 

NBR 

NWGR 

NMDPR 

NMDVR 

NMDTR 

NFNVR 

NFNTR 

NVPW 


Card  5: 

GC 
GH 


Gas  generator  diffuser  efficiency 

Compressor  efficiency 

Burner  efficiency 

Turbine  efficiency 

Gas  generator  nozzle  efficiency 

Fan  duct  diffuser  efficiency 

Fan  efficiency 

Fan  nozzle  efficiency 

Water  injection  pump  efficiency 


Ambient  density  (slug/fr  ) 

Ambient  pressure  (lb/ft  ) 

Ambient  temperature  (degrees  Rankine) 

Heating  value  of  fuel  (BXiJ/lbm) 

Reference  air  temperature  (492  deg,  R) 


Air  injection  mach  number  (at  fan  duct  entrance) 

Sea  water  temperature  (deg,  R) 

Burner  total  pressure  ratio 

Fan  duct  total  pressure  ratio  (dry  fan) 


Number  of  vessel  speeds 
Number  of  compressor  pressure  ratios 
Number  of  fan  pressure  ratios 
Number  of  turbine  inlet,  temperatures 
Number  of  bypass  ratios 

Number  of  water-to-gas  generator  air  ratios 
Number  of  mixing  duct  pressurt  rauos 
Number  of  mixing  duct  velocity  ratios 
Number  of  mixing  I  i  '  :.t-mperaf  ire  radios 
Number  of  fan  nozzle  velocity  ratios 
Number  of  fan  nozzlt  tempera t  :re  ratios 
Number  of  water  injection  velocities 


Cool  temperature  specific  hear  ratio 
Hot  temperature  specific  hear  ratio 
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Card  6: 
VO(I) 

Card  7: 

PTR21(J) 
Card  8: 

PTR76(M) 
Card  9: 

TT3(K) 
Card  10: 

BR(L) 
Card  11: 

WGR(N) 
Card  12: 

PSR87(IJ) 
Card  13: 

VR8WA(IK) 
Card  14: 

TSR8WA(IL) 
Card  15: 

VR9WA(IM) 
Card  16: 

TSR9WA(IN) 
Card  17: 

VPW(NN) 


Values  of  vessel  speed  (knots) 


Values  of  compressor  pressure  ratio 


Values  of  fan  pressure  ratio 


Values  of  turbine  inlet  temperature  (deg.  R) 


Values  of  bypass  ratio 


Values  of  water-to-gas  generator  air  ratio 


Values  of  mixing  duct  static  pressure  ratio 


Values  of  mixing  duct  velocity  ratio 


Values  of  mixing  duct  temperature,  ratio 


Values  of  fan  nozzle  velocity  ratio 


Values  of  fan  nozzle  temperature  ratio 


Values  of  water  injection  velocity  (knots) 


OUTPUT  PARAMETERS  NOT  PREVIOUSLY  DEFINED 

PWORK  Water  injection  pump  work  (ft- lb/slug) 

DHT43  Turbine  work  (ft-lb/slug) 

STHRUS  Specific  thrust  (lb  thrust/lbm/sec ) 

TR  Ratio  of  water-augmented- to-dry  specific  thrust 

ETAPRO  Propulsive  efficiency 

ERATIO  Ratio  of  water-augmented- to-dry  propulsive  efficiency 


79 


SFC 
VA7 
TS8 
VA8 
TS9 
VA9 


Specific  fuel  consumption  (lb  fuel/lb  thrust-hr) 

Air  velocity  at  fan  duct  entrance  (ft/sec) 

Air  static  temperature  at  mixing  duct  exit  (deg.  R) 

Air  velocity  at  mixing  duct  exit  (ft/sec) 

Air  static  temperature  at  fan  nozzle  exit  (deg.  R) 

Air  velocity  at  fan  nozzle  exit  (ft/sec) 
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APPENDIX   III 

DATA  REDUCTION  PROGRAM  NOMENCLATURE 

INPUT  PARAMETERS 

M_        Number  of  interpolation  pairs  for  subroutine  SPLIN1 
R_  _  (I)   Reynolds  number  values 

K_  _  (I)   Discharge  coefficient  values  corresponding  to  R_   (1) 
Subscripts  for  M_  _,  R__  _  (I),  K (I): 

FA    Flange  pressure  taps,  air 

DA    D  and  \   D  pressure  taps,  air 

FW    Flange  pressure  taps,  water 

DW    D  and  %  D  pressure  taps,  water 
Example:   RDA(I)   Reynolds  number  for  D  and  \   D  taps  for  air  flow 

M_        Number  of  interpolation  pairs  for  subroutine  SPLIN1 
X_  (I)    Lockhart-Martinelli  two-phase  flow  modulus  values,  ^C 
F_     (I)  Lockhart-Martinelli  correlation  parameter  values,  Q 
Subscripts  for  M_,  X_(I),  F_  _  _  (I): 

G     Lockhart-Martinelli  gas  correlation 
L     Lockhart-Martinelli  liquid  correlation 
TT    Turbulent- turbulent  flow 
TV    Turbulent-viscous  flow 
VT    Viscous- turbulent  flow 
W    Viscous-viscous  flow 
Example:   FGTV(I)   Correlation  parameter,  <Pq,    for  gas-phase 
turbulent-viscous  flow 

D1A  Inside  diameter  of  air  pipe  at  orifice  (in) 

D2A  Diameter  of  air  orifice  (in) 

D1W  Inside  diameter  of  water  pipe  (in) 

D2W  Diameter  of  water  orifice  (in) 

HMFA  Manometer  fluid  differential,  flange  taps,  air  fin  HO) 

HMDA  Manometer  fluid  differential,  D  and  \   D  taps,  air  (in  HO) 

HMFW  Manometer  fluid  differential,  flange  taps,  water  (in  Hg ) 

HMDW  Manometer  fluid  differential,  I)  and  \    D  taps,  water  (in  Hg) 
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T1FA      Air  flow  temperature  at  orifice  (deg.  F) 

T1FW      Water  flow  temperature  at  orifice  (deg.  F) 

TAFA      Air  temperature  at  air  orifice  manometer  (deg.  F) 

TAFW      Air  temperature  at  water  orifice  manometer  (deg.  F) 

2 
PI        Air  static  pressure  upstream  of  orifice  (lb/in  ) 

2 
PT        Static  pressure  at  test  section  inlet  (lb/in  ) 

TTFA  Air  temperature  at  tesn  section  outlet  (deg.  F) 

TTFW  Water  temperature  at  test  section  outlet  (deg.  F) 

DT  Inside  diameter  of  test  section  (in) 

LT  Length  of  test  section  (ft) 

OUTPUT  PARAMETERS  NOT  PREVIOUSLY  DEFINED 

BETAA  Ratio  of  orifice-to-pipe  inside  diameter,  air  pipe 

BETAW  Ratio  of  orifice-to-pipe  inside  diameter,  water  pipe 

RFAI  Air  Reynolds  number  based  on  flange  taps 

RDAI  Air  Reynolds  number  based  on  D  and  \   D  taps 

RFWI  Water  Reynolds  number  based  on  flange  taps 

RDWI  Water  Reynolds  number  based  on  D  and  \   D  taps 

KFA1  Air  discharge  coefficient,  based  on  flange  taps 

KDAI  Air  discharge  coefficient  based  on  D  and  %  D  taps 

KFWI  Water  discharge  coefficient  based  on  flange  taps 

KDWI  Water  discharge  coefficient  based  on  D  and  \   D  taps 

FRFA  Air  mass  flow  rate  based  on  flange  taps  (lbm/sec) 

FRDA  Air  mass  flow  rate  based  on  D  and  \   D  taps  (lbm/sec) 

FRFW      Water  mass  flow  rate  based  on  flange  taps  (lbm/sec) 

i 

FRDW      Water  mass  flow  rate  based  on  D  and  \   D  taps  (lbm/sec) 

D         Inside  diameter  of  test  section  (ft) 

3 
RHOWT     Water  density  in  test  section  (lbm/ft  ) 

3 

RHOAT     Air  density  in  test  section  (lbm/ft  ) 

RLM       Liquid  Reynolds  number  (Lockhart-Mar tinell i) 
RGM       Gas  Reynolds  number  ( Lockhar  t.-Mart  inelli  ) 
VL        Liquid  velocity  (ft/sec) 

VG        Gas  velocity  (ft/sec) 

2 
PDROPL    Fictitious  liquid,  pressure  drop  (lb/ft  /10ft) 

2 
PDROPG    Fictitious  gas  pressure  drop  (lb/ft  /10ft) 

TPPDi     Lockhar t -Mar tine Hi  two-phase  pressure  drop  prediction 
(lb/in  /10ft) 
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RLC  Liquid  Reynolds  number  ( Che  nerve  th- Mart  in) 

RGC  Gas  Reynolds  number  (Chenoweth-Mar tin) 

PDLC  Dimensionless  group s  ~ty       (all- liquid  flow) 

PDGC  Dimensionless  group,  al*     (all-gas  flow) 

CORPAR  Chenoweth-Mar tin  correlation  parameter 

VFR  Volumetric  flow  ratio  (liquid-to-gas) 

2 

PDRLC  Fictitious  all-liquid  pressure  drop  (lb/ft  /10ft) 

LVOLFR  Liquid  volume  fraction 


DATA  INPUT  CARD  ORDER 


Card  1 

:    MFA,  MDA,  1 

MFW, 

MDW 

Card  2 

:    RFA(I), 

I 

■  1, 

MFA 

Card  3 

:    KFA(I), 

= 

1,  MFA 

Card  4 

:    RDA(I), 

I 

=  1, 

MDA 

Card  5 

KDA(I), 

I 

=  1, 

MDA 

Card  6 

'    RFW(I), 

I 

■  1", 

MFW 

Card  7 

KFW(I), 

I 

=  1, 

MFW 

Card  8 

RDW(I), 

I 

■  K 

MDW 

Card  9 

KDW(I), 

1 

=  i, 

MDW 

Card  1( 

):   D1A,  D2A, 

D1W, 

D2W 

Card  1] 

_:   HMFA,  HMDA 

,  T1FA,  TAFA,  PI 

Card  U 

!:   HMFW,  HMDW 

,  T1FW,  TAFW 

Card  i: 

5:   PT,  TTFW, 

FTFA 

Card  U 

v:      DT,  LT 

Card  If 

i:   MG,  ML 

Card  16 

.:   XG(I),  : 

[  = 

1,  MG 

Card  17 

:  XL(I),  : 

[  = 

1,  Ml. 

Card  U 

:   FGTT(I) 

I 

=  1, 

MG 

Card  IS 

>:   FGTV(I)_ 

I 

«  1, 

MG 

Card  2C 
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=  1 

MG 

Card  21 
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■  I, 

MG 

Card  21 

:   FLTT(I), 

I 

=  1: 

ML 

Card  2: 
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■  1, 

ML 

Card.  24 

r.      FLVT(I), 

I 

-  1  = 

ML 

Card  2f 

:   FLW(I), 

I 

-  1. 

Ml, 
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several  variables,  including  water-to-gas  generator  air  ratio,  water  injection  veloc- 
ity, fan  duct  pressure  loss,  and  fan  duct  thermal  and  dynamic  nonequilibrium,  upon 
thrust  and  propulsive  efficiency.   In  addition,  the  design  parameters  of  fan  pressure 
ratio  and  fan  bypass  ratio  are  ^examined  for  their  optimum  values,  and  optimum  operat- 
ing combinations  of  water-to-gas  ratio  and  water  injection  velocity  are  determined. 

A  test  apparatus  is  developed  for  the  direct  measurement  of  wall  friction  force 
in  two-phase  flows.  A  computer  program  is  presented  to  reduce  experimental  data  and 
compare  with  pressure  drop  predicted  by  two  empirical  correlations. 
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